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RESEARCH  AT  NRL  ON  THEORETICAL  AND  NUMERICAL 
SIMULATION  STUDIES  OF  IONOSPHERIC  IRREGULARITIES 
!•  Introduction 

The  general  problem  of  Ionospheric  modeling,  important  for  pre- 
dictive capabilities  in  communication  and  radar  systems,  can  be  broken 
up  into  three  categories:  (1)  data  collection;  (2)  model  development; 
and  (5)  systems  applications.  The  Plasma  Dynamics  Branch,  the  theo- 
retical and  computational  arm  of  the  Plasma  Physics  Division,  has  been 
engaged  in  the  second  category  of  dynamic  ionospheric  modeling  for 
about  two  years.  This  research  is  being  carried  out  under  the  auspices 
of  ONR,  NAVAIR,  and  DNA.  The  general  modeling  categories  can  be 
broken  up  into  the  following  subsections:  (i)  chemistry;  (ii)  neutral 
atmosphere;  (lii)  laminar  ionosphere;  and  (iv)  turbulent  or  Irregular 
ionosphere.  The  work  which  will  be  reported  on  in  this  paper  is  an 
overview  of  the  studies  carried  out  in  area  (iv)  which  we  call  iono- 
spheric plasma  irregularities.  Study  of  the  turbulent  ionosphere 
requires  a knowledge  of  the  laminar  ionosphere  which  in  turn  requires 
a knowledge  of  the  chemistry  and  the  neutral  atmosphere.  Thus,  for  a 
time  predictive  capability  areas  (1)  - (iv)  must  be  studied.  A program 
addressing  all  of  these  areas  is  being  developed.  The  present  report 
will,  however,  address  only  the  problem  of  ionospheric  irregularities. 

Initially  the  work  involved  the  development  of  the  capability  to 
treat,  computationally  and  theoretically,  the  non-linear  dynamics  of 
ionospheric  plasmas.  The  capability,  once  developed,  was  then  applied 
to  several  outstanding  problems.  Three  of  the  problems  are:  (A)  Type 

II  equatorial  electro jet  instabilities;  (B)  Striations  in  plasma  clouds 
and  image  formation  in  plasma  cloud- ionosphere  coupled  systems;  and  (C) 
Ionospheric  spread  F.  Radar  backseat ter^  observations  of  three  meter 
Manuscript  submitted  September  27,  1974. 
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irregularities,  i.e.,  electron  density  fluctuations,  in  the  geomag- 
netic equatorial  E region  electrojet  have  led  us  to  study  process  A. 

2 5 

Barium  clouds  * released  by  rockets  in  the  upper  ionosphere  and  high 
altitude  nuclear  releases  both  of  which  exhibit  field  aligned  striations, 

i.e.,  electron  and  ion  density  fluctuations,  led  us  to  the  study  of 

k 

process  B.  Scintillation  of  radio  signals  from  polar  orbiting  and 

56  789 

equatorial  satellites  ’ and  ionosonde  data  * have  led  us  to  study 
process  C (also  due  to  electron  density  fluctuations).  All  three 
processes  can  cause  scintillation  or  enhanced  backscatter.  These  three 
processes  are  caused  by  plasma  instabilities  which  depend  on  the  fact 
that  electrons  and  ions  collide  with  neutral  particles.  In  particular, 
the  instabilities  causing  the  first  two  processes  (A  and  B)  arise  from 
having  an  electron  density  gradient  and  electron  and  ion-neutral  drag 
effects  transverse  to  the  ambient  magnetic  field.  This  instability  is 
called  the  E X B gradient  drift  instability. Process  C, 
in  the  case  of  midlatitude  Spread  F,  is  believed  to  depend  upon  the 
fact  that  as  the  plasma  in  a flux  tube  is  moved  up  or  down  in  altitude 
it  samples  regions  of  different  Pedersen  conductivity^^.  The  inadequacy 
of  linear  theory  to  describe  these  phenomena  has  led  us  to  study  by 
theoretical  and  numerical  simulation  techniques  the  nonlinear  develop- 
ment of  these  three  processes  (plasma  instabilities). 
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II.  Importance  of  Studies 


The  basic  importance  of  our  studies  on  ionospheric  irregularities 
can  be  catalogued  as  follows:  (1)  the  instabilities  degrade  comnunica- 

tions;  (2)  the  instabilities  degrade  radar  performance;  and  (3)  the 
studies  will  enable  us  to  obtain  a predictive  capability  with  respect 
to  (1)  and  (2),  i.e.,  where  to  expect  these  Irregular  structixres  and 
what  their  properties  will  be.  At  this  juncture  let  us  elaborate  on 
(1). 

The  fact  that  time-varying  irregular  structure  in  the  ionosphere 
(known  as  ionospheric  scintillation)  could  affect  radio  signals  emit- 
ted  by  stars  has  been  known  for  some  time  ’ . However,  only  recently, 

and  then  with  some  surprise,  has  it  been  found  that  the  ionosphere  can 
seriously  perturb  signals  used  for  satellite  communications^'^.  The 
effect  of  Irregular  structure,  i.e.,  irregularities  of  the  refractive 
index,  in  the  ionosphere  is  to  superpose  a random  fluctuation  of  signal 
amplitude,  phase  and  polarization.  These  deleterious  effects  may  cause 
an  intolerable  number  of  errors  in  the  received  signal  and  produce 
unusable  information,  Utlaut  has  recently  pointed  out  that  (a)  some 
of  the  orbiting  meteorological  satellites  transmitting  video  vreather 
data  have  experienced  conditions  during  which  the  video  data  were  near- 
ly obscured  as  a result  of  scintillation;  (b)  a navigation  satellite 
system  has  experienced  difficulties  inserting  ephemeris  data  in  the 
spacecraft  at  times  because  of  scintillation;  and  (c)  signals  from 
geostationary  satellitles  may  also  suffer  deleterious  effects,  in  regions 
other  than  the  equatorial  region,  when  transmitting  to  higher  latitudes 
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since  the  polar  cap  scintillating  region  is  tipped  with  respect  to  the 
geographic  pole  (e.g,,  aircraft  routes  near  midlatitude  over  the  North 
Atlantic  could  experience  difficulties  communicating  by  way  of  geosta- 
tionary satellites).  Naturally  occurring  scintillation  is  strongest  in 

the  high  latitude  (auroral  and  polar  regions)  and  equatorial  regions^"^’ 

18  IQ  20  21 

, although  midlatitude  scintillation  also  occurs  * . However, 

these  notions  refer  to  natural  scintillations.  U.  S.  high  altitude 

nuclear  releases  (e.g.  Starfish)  have  produced  world  wide  spread  F with 

22  25  24 

attendant  scintillation  phenomena  * * . Also,  Russian  releases  at 

Novaya  Zemlya,  U.S.S.R.  in  October  I96I  produced  E and  F region  irregu- 

25 

larltles  in  the  Scandinavian  area  . 

Although  much  of  the  scintillation  observations  have  been  restrict- 

17 

ed  to  VHF  , recent  results  of  equatorial  scintillation  in  the  1-6  GHz 

^ X *7 

regime^’  * show  that  one  cannot  merely  go  to  higher  frequencies  to 
escape  scintillation  effects.  Much  of  the  scintillation  observations 
come  from  ground  based  receivers  (mostly  in  VHF  and  UHF  regime,  with  a 
small  amount  of  SHF  data)  although  limited  in  situ  measurements  are 
also  available.  These  observations  have  provided  valuable  information 
about  scintillation  phenomena;  however,  the  effects  of  variation  of 
signal  parameters  important  for  system  design  or  systems  performance 
prediction  is  almost  nonexistent.  More  fundamental  is  the  lack  of  under- 
standing of  the  basic  physical  processes  involved  in  the  formation  of 
ionospheric  irregularities  which  cause  scintillation.  The  limited 
amount  of  observational  data  coupled  with  limited  theoretical  under- 
standing make  it  impossible  to  now  predict  how  communications  systems 
will  perform  in  natural  or  artificial  scintillation  environments.  For 
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the  high  data  rate,  broadband  transmitters  now  planned  for  Navy,  other 
DOD  and  commercial  communication  satellites  it  is  imperative  to  under** 
stand  scintillation  phenomana,  for  severe  scintillations  over  several 
• seconds  can  mean  loss  of  Information  (message) • 

Ionospheric  physicists,  systems  designers,  modelers  and  wave  propaga- 

I7  26 

tionists  have  noted  * the  need  for  theoretical  studies  to  determine 
the  underlying  mechanisms  responsible  for  ionospheric  scintillation  so 
that  a true  predictive  capability  can  be  achieved.  The  inadequacy  of 

26 

empirical  modeling  of  ionospheric  scintillation,  the  recent  in  situ 

21 

measurements  of  ionospheric  scintillation  , and  advances  in  theoretical 
and  computational  plasma  physics^^*^^"^^  make  this  a propitious  time  for 
further  theoretical  and  numerical  studies  of  scintillation  phenomena. 

w 

Ihe  theoretical  and  computational  programs  at  NRL  are  at  a state,  where 
^ we  can  generate  detailed  Information  on  the  causes  and  structure  of 

ionospheric  irregularities.  The  results  of  this  work  can  provide  a 
basis  for  interpretation  of  experimental  data;  suggest  new  experiments; 
aid  in  the  design  of  new  communication  systems  and  provide  a framework 
within  which  to  improve  the  performance  of  existing  systems. 


III.  Simplified  Scintillation  Geometry 

Before  proceeding  with  a description  of  examples  of  the  work  we 
are  doing  on  ionospheric  irregularities,  we  will  give  a brief  discussion, 
albeit  simple,  on  what  is  involved  in  scintillation  phenomena.  The 
interested  reader  can  obtain  a more  complete  understanding  of  propagation 
through  ionospheric  irregularities  by  consulting  the  appropriate 
references^*^*^^”^^ , 
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Figure  1 shows  a simplified  view  of  what  happens  in  scintillation 
phenomena.  We  consider  a plane  wave  (signal)  propagating  in  the  * 
direction  where  it  is  emitted  by  a transmitter  (source)  and  enters  a 
layer  where  there  are  irregularities,  i.e.,  density  fluctuations  • 

(shown  by  ellipses  in  the  x-y  diffracting  plane).  On  emergence  from 
the  irregularities  layer,  properties  of  the  wave  front  are  related  to 
properties  of  the  irregularities.  In  most  cases  of  Interest  absorption 
in  the  layer  is  negligible.  Also,  the  irregularity  layer  can  usually 
be  replaced  by  an  equivalent  thin  screen  of  negligible  thickness. 

Field  fluctuations  near  the  screen  are  of  phase  only;  however,  anq>litude 
scintillations  develop  as  the  wave  travels  away  from  the  screen.  For 
these  conditions  the  relationship  between  the  intensity  scintillation 
spectrum,  (k),  and  the  irregularities  density  fluctuation  spectrum, 

Pjj  (k)  is^^ 

Pg  (k)  = 8ttL  (r^X)^  Pjj  (k)  sin^  (k/k^)^  (1) 

where  L is  the  equivalent  layer  thickness,  r is  the  classical  radius 

e 

of  the  electron  (r^  = e^/mc^),  X is  the  radio  wavelength  of  the  signal, 

k is  the  wavenumber  of  the  irregularity  of  interest  (k^  = k ^ + k ^) , 

X y ’ 

k^  is  the  Fresnel  wavenumber  (k^^  = 4tt/Xz,  where  z is  the  mean  height 
of  the  scattering  layer) . Equation  1 gives  the  expected  two-dimensional 

diffraction  pattern  at  the  receiving  plane.  Indeed  the  integral, 

00 

Pg  (k)  k d k is  basically  the  scintillation  index  used  by 
ionospheric  workers.  The  important  point  is  that  there  exists  a simple 
relationship  between  Pg  (k)  and  the  power  spectrum  P^^  (k)  of  the 
irregularities  (even  for  thick  screens;  see  reference  20).  P„  (k) 
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determines  the  diffraction  pattern  and  (k)  determines  (k) . Thus 
vhat  we  wish  to  find  out  from  the  theoretical  and  numerical  simulation 
studies  Is  the  power  spectrum,  P^  (k),  of  the  density  fluctuation^  among 
other  things.  To  date,  P^  (k)  is  assumed  by  Ionospheric  workers  (based 
on  at  best  limited  data  on  scintillation  phenomena) . To  obtain  P^  (k) 
from  theoretical  and  ntmerical  simulation  analysis  will  give  a true 
predictive  capability. 


IV.  Particular  Examples  of  Irregularities  Studied 

Before  proceeding  to  discuss  our  work  on  the  type  II  electrojet 
Instability  and  strlatlons,  as  exanq>les.  In  sections  A and  B below 
(work  on  midlatitude  spread  F is  In  the  beginning  stages),  we  will 
briefly  outline  the  methods  of  solution  which  are  involved  and  basically 
common  to  all  three  processes  A,  B,  and  C mentioned  In  Section  I.  The 
methods  involved  are  (1)  large  scale  numerical  simulation  of  the  appro- 
priate system  of  coupled  nonlinear  differential  equations  which  are 
derived  theoretically  and  thought  to  represent  the  process  being 
studied;  (2)  the  geometry  is  basically  two-dimensional,  i.e.,  the 
dimension  parallel  to  the  ambient  magnetic  field,  j^.  Is  Integrated  out; 
(5)  there  are  continuity  eqxiatlons  for  the  density  and/or  cross  field 
(Pedersen)  conductivity  of  the  appropriate  electron  and  Ionic  species; 
(4)  momentum  equations  for  the  appropriate  species;  and  (^)  a potential- 
like equation  arising  from  the  fact  that  the  divergence  of  the  current 
is  zero  (v  • ^ • 
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The  two  dimensional  geometry  is  a fairly  good  representation  of 
the  processes  being  studied.  Full  three  dimensional  studies  would  tax 
the  capability  of  present  day  computers.  The  momentum  equation  in  (h) 
usually  reduces  to  an  algebraic  equation  In  which  the  velocity  of  the 
particular  species  Is  related  to  the  transverse  electric  fields  (and/or 
neutral  winds)  and  ambient  magnetic  field  by  means  of  cross  field 
mobilities.  Assuming  the  electric  fields  derivable  from  a potential 
(the  electrostatic  approximation  valid  for  the  Ionosphere)  results  In  a 
potentlal-llke  equation  when  the  velocities  are  Inserted  Into  the  current 
condition,  y • J = 0.  All  these  Ideas  will  become  clearer  In  sections  A 
and  B below, 

A,  Equatorial  Electro  let  - Type  II  Instability 
To  serve  as  a background  for  this  section.  Figure  2 shows  a model 
of  the  daytime  electro jet  (see  references  42  and  45)  and  how  the  radar 
views  It,  The  electro jet  Is  a highly  conducting  layer  usually  between 
100  and  110  km  located  at  the  geomagnetic  equator  (within  ± 2°) . The 
eastward  field,  , results  from  atmospheric  tidal  motion  and  the 
vertical  electric  field,  Is  the  resulting  vertical  polarization 
field  which  suppresses  vertical  currents,  (This  condition  gives 
Ep  = (ajj/op)  Ejj,  where  and  a-p  are  the  Hall  and  Pedersen  conduc- 
tivities respectively  and  a^/cp  20-50),  It  Is  this  vertical  field 
which  drives  the  electro jet  by  Hall  currents.  Since  the  electron- neutral 
collision  frequency  Is  much  less  than  the  electron  cyclotron  frequency 
(Vq  <<  Qg)  and  the  Ion-neutral  collision  frequency  Is  much  greater  than 
the  Ion  cyclotron  frequency  (v^>>  Ions  are  essentially  at  rest  with 
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respect  to  the  neutrals),  the  resulting  electron-ion  drift  velocity  is 

cEp 

Vq  « g — (cgs  units  will  be  used  throughout)  and  in  the  westward  direc- 
o 

tion*  This  is  essentially  the  electron  drift  so  that  the  current  is  in 
the  eastward  direction.  Also  there  is  a vertical  electron  density 
gradient  with  a scale  height  (N^/dN^/dz)  ~ 6-10  km.  The  reason  why  the 
electrojet  is  localized  in  altitude  is  that  above  this  region  while  N 

o 

is  increasing  the  collision  frequency  is  decreasing  and  so  the  conduc- 
tivity is  low.  However,  below  this  region  while  the  collision  frequency 
is  increasing  the  electron  density  is  decreasing  and  so  again  the  conduc- 
tivity is  low.  The  reason  why  the  electrojet  is  located  at  the  geo- 
magnetic equator  is  deduced  by  considering  the  concept  of  a stratified 
ionosphere  and  the  concept  of  layered  conductivity^^.  Basically  the 
requirement  that  there  be  no  vertical  current  reduces  the  5X5  conduc- 


tivity tensor  to  a 2 X 2 tensor  representing  a layered  conductivity  with 
horizontal  flow.  Into  such  an  expression  comes  the  dip  angle  (angle 
between  the  magnetic  field  and  the  horizontal  direction).  It  can  be 
shown  that  when  the  dip  angle  is  zero  the  cross  field  conductivity  be- 
comes a maximum.  Referring  to  Fig,  2 the  radar  observes  the  electrojet 
by  backscatter  so  that  a wave  launched  with  wavevector  ^ sees  a return 
with  -1^  and  the  relation  between  the  radar  wavelength  and  the  irregu- 


wavelength  is  2 |^j  = |^j|  where  k =s  27t/\  and  \ is  the  wavelength. 
(So  that  the  radar  wavelength  is  twice  the  wavelength  of  the  irregularity 
that  it  samples).  The  radar  at  Jicamarca,  Peru  where  the  observations 
are  made  operates  in  the  I5  MHz-150  MHz  range  (Irregularity  wavelengths 
10-1  meter) . 
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There  are  two  types  of  instabilities,  producing  density  fluctua- 
tions, associated  with  the  electrojet  configuration  shown  in  Fig.  2. 

In  this  paper  we  will  be  addressing  the  type  II  instability.  This  is 
excited  when  the  electron  current  across  ^ is  less  than  the  ion  sound 
speed  (i.e.,  < ion  sound  speed).  This  instability  has  no  velocity 

threshold  and  travels  at  the  drift  velocity  of  the  electrons.  The 

observed  backscatter  is  generally  believed  to  be  due  to  this  gradient 
13  1I4. 

drift  instability  ^ '.  However,  we  will  see  shortly  that  linear 

instability  theory  doesn't  bear  this  out  for  the  radar  wavelengths  ob- 
served, The  basic  equations  governing  the  dynamics  of  the  electrojet 


are 


52,55 


^ + V . (N^V  =0 


(2) 


U 


~ (E  + V X B ) + -r — y N + j/  V = 0 
m ~ “e  ^ N e e ~e 


(5) 


^ + V . (N . V ) = 0 

at  ^ ^ 


(4) 


=-57’  + S; £ - •'I 


(5) 


N = N = N 
e i 


(6) 


where  the  subscripts  e and  i refer  to  electrons  and  ions,  U and  w 

* e *^e 

are  the  electron  thermal  speed  and  electron-neutral  collision  frequency. 
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respectively,  E = - y cd  is  the  electric  field  (electrostatic  assumption), 
is  the  earth's  magnetic  field  (assumed  constant  and  in  the  x direc- 
tion), N is  density  and  V is  bulk  velocity.  The  inclusion  of  the  left 
hand  side  eqn.  (5)>  i#e.,  ion  inertia,  is  not  really  necessary  for  type 

II  studies.  However,  it  is  important  for  type  I studies  (V  > ion 

o 

sound  speed)  and  we  retain  it  for  generality.  Equation  (3)  contains 

magnetic  terms  (j;  « n ) showing  the  electrons  are  tied  to  the  field 
e e 

lines,  while  eqn  (5)  exhibits  no  such  terms  » fj^)  showing  the  ions 
to  be  tied  to  the  neutrals. 

A linear  stability  analysis  of  eqns.  (2)-(6)  with  perturbations 

assumed  to  be  of  the  form  exp  i (k^  y - ujt),  where  cu  = + i Y (so 

that  positive  y represents  the  growth  rate),  has  been  performed.  For 

electrojet  parameter 0^=  5 X 10®  sec”^,  = 4 X 10“*  sec”^, 

V.  = 2.5  X 10®  sec  C ^ = 10®  m^/sec^,  (where  C_^  = U.^  + U ^ m /m, ) 

1 b b 1 e e 1 

V = 100  m/sec  (V  < C^,  ion  sound)  and  (N  /dN  /dz)  = 6 X lO^m.  unstable 

O Ob  oo 

gradient  drift  waves  result  in  the  y direction  with  a minimum  wave- 
length of  50  meters  (see  Fig.  5).  However,  radar  backscatter  results^ 
show  unstable  waves,  i.e.,  enhanced  density  fluctuations  3 meters. 

In  order  to  resolve  this  dilemma,  we  have  studied  the  nonlinear  evolution 

of  the  electrojet  equations  (2)  - (6),  in  the  y-z  plane,  by  numerical 
33 

simulation  . 

Before  proceeding  with  the  numerical  simulation,  eqns.  (2)  - (6) 
are  recast  using 

v.j_  = e V . [ N (V^  - V^)J  (7) 

As  stated  previously  eqn  (7)  results  in  a Poisson-like  equation  by 


11 


algebraically  solving  eqn.  (3)  for  and  using  an  effective  potential, 

«)  = — 0 - U ^ In  N (8) 

m ® 

e 

The  resulting  equation  for  d>  is 

o r -e  1 1 

V^<I.+  |vlnN-  — XvlnNj.v*=(v^+  — )^  y.(N  V^)  (9) 

where  O = e B /m  c. 

•^e  — o e 

An  initial  density  perturbation^^ of  the  form  6 N/N  = 0.05 

o 

sin  ky  y - 0.01  sin  z,  where  k^  = k^  = 2^7/72  meters  is  put  into  the 
system  at  t = 0 (see  Fig.  4).  We  are  essentially  looking  at  a slice  of 
the  electro jet  in  the  y-z  plane  which  is  75  x 75  meters.  We  then 
numerically  integrate  eqns.  (4),  (5)  and  (9)  and  carry  it  forward  in 
time  (see  reference  55  for  a more  complete  discussion  of  the  procedure 
used)  and  what  results  is  a quasi-final  state  of  turbulence  as  depicted 
in  Fig.  5. 

The  density  contours  depicted  in  Fig.  5 Is  the  type  of  turbulent 
layer  that  a signal  traversing  the  E region  electrojet  would  see.  The 
relevant  Information  from  such  a figure  is  obtained  by  performing  a 
Fourier  analysis  of  the  ( 6 N/N)^  and  this  power  spectrum,  i.e.  the 
(k)  we  spoke  about  in  section  III,  is  shown  in  Fig.  6.  We  note 
that  this  square  of  the  density  fluctuation  spectral  amplitude  follows 
a power  law  for  the  shorter  wavelengths  which  agrees  well  with  rocket 
observations  through  the  equatorial  electro jet^®**^^ . It  is  of  interest 
to  note  that  Spread  F in  sltu^^  and  ground  observed^^’^^  scintillation 
measurements  show  similar  power  law  dependence  of  the  spectrum. 
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Type  II  electro jet  results  can  be  summarized  as  follows:  (1)  long 
wavelength  (>  50  meters)  horizontally  propagating  gradient  drift  in- 
stabilities generate  short  wavelength  (~5  meters)  vertically  propa- 
gating instabilities;  (2)  the  quasi-final  state  is  a highly  turbulent 
two  dimensional  state  with  upward  and  downward  moving  irregularities 
with  speeds  of  the  order  of  the  horizontal  drift  speeds;  and  (5)  the 
power  spectrxim  for  density  fluctuations  in  the  short  wavelength  irre- 
gularitles  goes  like  k Result  (1)  agrees  with  the  high  resolution 

radar  backscatter  data^  and  our  basic  theoretical  ideas^^*^^.  It  basi- 
cally shows  a cascading  of  energy  from  longer  linearly  unstable  wave- 
lengths to  shorter,  linearly  stable  wavelengths.  Result  (2)  also  agrees 
well  with  the  radar  backscatter  data^.  Result  (5)  agrees  well  with  the 

i<.8  4q  so 

rocket  observations  ’ ^ and  dimensional  arguments^  . 

This  is  the  first  numerical  simulation  ever  done  of  the  type  II 
electrojet  phenomena  and  has  aided  considerably  in  the  understanding  of 
such  phenomena.  Further  simulations  are  on  going  to  obtain  finer 
resolution  for  the  shorter  wavelengths  (see  Fig.  6 and  comments  therein) 
and  also  Increasing  so  that  the  intermediate  regime  between  type  II 
and  type  I can  be  studied.  The  capability  which  has  been  developed  for 
this  problem  can  also  be  used  to  study  the  stability  of  the  auroral 
electro jets. 
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B,  Strlations  and  Cloud  Morphology  Associated  With  Plasma 


Clouds  Coupled  to  the  Background  Ionosphere 

We  will  now  discuss  our  two-layered  model  for  striation  and  image 

cloud  development  when  plasma  clouds  exist  in  and  are  coupled  to  the 

background  ionosphere.  Striations  are  magnetic  field-aligned  irregu- 

2 5 

larities  such  as  are  seen  in  barium  plasma  clouds  ^ released  in  the 
ionosphere  and  in  the  late  time  development  of  high  altitude  nuclear 
explosions.  These  phenomena  are  often  diagnosed  by  optical  techniques. 
If,  however,  these  striations  induced  irregularities  in  other  regions  of 
the  ionosphere  such  irregularities  would  be  invisible  to  the  naked  eye, 
but  not  to  radar  or  communication  signals  that  pass  through  them.  In- 
deed, these  image  effects,  as  we  shall  call  them,  could  have  a profound 
effect  on  such  signals.  Although  we  have  spoken  of  artificially  released 
plasma  clouds  inducing  these  image  effects,  the  basic  ideas  are  equally 
applicable  to  naturally  occurring  plasma  blobs  or  irregularities  that 
appear  in  the  ionosphere  (here  the  clouds  themselves  are  only  visible 
to  radio  signals'.).  Figure  7 shows  a simplified  view  of  our  basic  two- 
level  system  for  treating  plasma  clouds  coupled  to  the  background  iono- 
sphere. Layer  1 represents  the  level  at  which  the  basic  plasma  cloud 
exists  and  layer  2 is  the  level  at  which  images  will  occur  (in  actuality 
images  will  extend  continuously  below  the  cloud  to  some  level  depending 
on  background  ionospheric  ion-neutral  collision  frequencies  and  ratio 
of  cloud  to  background  ionosphere  Pedersen  conductivity).  Each  of  the 
two  levels  is  characterized  by  its  own  zeroth  order  parameters  such  as 
ion- neutral  collision  frequencies,  cross  magnetic  field  conductivities. 
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etc.  (see  for  example  references  50  and  3I) . The  two  regions  are 
electrostatically  coupled  because  of  the  high  conductivity  along  the 
magnetic  field  lines  which  couple  the  two  layers.  Tjrpically  we  consider 
cases  where  the  plasma  cloud  is  in  the  F region  and  coupled  to  the  E 
region,  i.e.,  the  cloud  is  around  200-250  km  and  would  produce  images 
down  to  around  I50  km.  Figure  7 gives  a pictorial  view  of  the  plasma 
cloud  in  layer  1 with  say  a gausslan  density  dependence  and  the  attend- 
ant image  clouds  induced  by  it  in  layer  2.  The  plus  and  minus  sign  in 
layer  2 refer  to  enhancement  of  positive  ions  and  depletion  of  positive 
ions  from  the  ambient  values,  respectively.  The  geometry  of  the 
applied  electric,  E^,  and  magnetic,  B^,  fields  are  shown.  Perturbations 
(striations)  produced  in  the  plasma  cloud  will  also  be  exhibited  in  the 
image  clouds^^’^^.  In  the  model  presented  in  Fig.  7,  ion  flow  is 
allowed  in  each  layer  and  essentially  only  electron  flow  between  layers 
(direction  depending  on  direction  of  E ^).  The  major  importance  is  the 
coupling  of  the  plasma  cloud  to  other  regions  of  the  ionosphere 
(imaging).  Indeed,  in  certain  cases,  one  can  get  coupling  to  the  con- 
jugate ionosphere. 

The  basic  starting  point  of  our  nonlinear  studies  of  striation 
phenomena  is  essentially  the  set  of  eqns.  (2)  - (6).  These  equations 
neglect  neutral  wind  effects  and  we  will  do  so  here,  initially.  It  is 
generally  agreed  upon  that  striatiore  in  plasma  clouds  are  caused  by 
the  E^  X gradient  drift  instability^^ and  so  the  left  hand  side 
of  eqn,  (5)  can  be  neglected.  However,  because  we  will  be  dealing  with 
altitude  regimes  where  the  ion-neutral  collision  frequency,  is  not 
much  greater  than  the  ion  cyclotron  frequency,  a magnetic  V X B force 
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term  must  be  included  in  eqn,  (5)  which  is  of  the  same  form  as  that 
which  appears  in  eqn.  (5) . For  example  the  perpendicular  ion  velocity 
can  then  be  solved  for  algebraically  in  terms  of  E and  B from  eqn  (5) 
and  we  obtain 


V - c \ Ti- V!i! 

' ®o  1..  - ^ ^ "i  -I 

^ (10) 

B = i«,^  L » »i  J 

O 1 

where  h Cl^/u^9  is  the  ion  cyclotron  frequency,  is  the  ion- 
neutral  collision  frequency,  and  j.  denotes  perpendicular  to  the  magnetic 
field  (which  in  this  study  is  in  the  z direction) . The  multiplicative 
factor  involving  in  eqn.  (10)  will  change  depending  on  what  altitude 
we  are  considering,  i.e.,  it  will  be  different  for  layers  1 and  2,  In 
this  study  we  point  out  that  there  are  basically  two  sets  of  ion 
equations  one  for  the  plasma  cloud  (which  we  refer  to  as  barium)  and 
one  for  the  background  ionospheric  ions  in  layer  2 (see  reference  50 
for  more  details).  The  quasi-neutrality  condition  eqn  (6)  becomes 


(11) 

where  i refers  to  background  ionospheric  ion  and  b to  barium  cloud  ion. 
We  take  the  electric  field  in  eqn,  (10)  to  be  jE  = E ^ ” V co»  where  E ^ 
is  ambient  field.  For  this  section  we  use  the  effective  potential  in 
the  form 


m U *= 

Y = 0 - e e 


In  N 


(12) 


and  is  constant  along  field  lines.  We  use  the  continuity  eqn  (4)  for 

barium  ions  and  integrate  over  the  entire  thickness  of  the  barium  cloud 
along  the  magnetic  field  line  and  obtain 
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b 


CVYXB^ 


= 


(lo-’f) 


vS 


(13) 


+ E.  Y + 2W  ^ p 

e e 


where  we  are  in  a frame  moving  with  c E X B / have  set  m U ^ = 

— O “O  r»  » A A 


kT 


o 

-1 


e e 


, have  kept  terms  to  first  order  in  where  refers  to 


barium  quantities  (this  is  possible  because  for  our  F region  cloud 
» 1) , is  now  only  a function  of  X,  y as  is  ^ which  is 


V,. 


T = ®£  N,  ^ 
P B 


_b 


(U) 


the  height-integrated  (field  line  integrated)  Pedersen  conductivity 
of  the  barium  cloud  (V|j/f)^  is  considered  to  be  constant  in  layer  1)  . 
Taking  the  limit  ->  <»  we  obtain  for  the  background  ionosphere  ions 


aN,  c vyxb 
1.  _ o 


at 


vN  ^ = - RHS  eqn.  (13) 


B 


(15) 


Equation  (7)  yields  after  intergration  over  the  z direction 
-(Sp%7'p'’)  v=  f +^E^- 4.Sp’’) 

where  the  superscript  i on  S refers  to  the  height  integrated  con- 
ductivity for  the  background  ionosphere  ions  and  have  used  the  fact 


(16) 


IT 


the  Hall  conductivity  of  the  ionospheric  ions 


H 


ecN^  = - I'i  E 


B 


1 + 

V 


(17) 


k) 


(note:  is  assumed  constant  in  layer  2;  also  for  a more  complete 

description  of  the  derivation  see  reference  50) • The  set  of  equations 

(15) ,  (15)  and  (16)  describe  the  two-level  system  depicted  in  Fig.  7 
and  are  two  dimensional  equations  (x,  y plane) , 

The  general  set  of  equations  describing  the  two  level  system  are 

not  exactly  transparent.  Consequently,  in  order  to  gain  a better 

understanding  of  the  physics  involved  in  the  coupled  cloud-ionosphere 

system,  we  have  broken  up  our  studies  into  small  and  large  clouds.  Here 

the  word  small  or  large  is  in  the  electrical  sense,  i.e.,  whether  the 

ratio  of  height-integrated  Pedersen  conductivity  of  the  cloud  to  that 

of  the  background  ionosphere  is  small  or  large.  The  small  cloud 
/ b i V 

approximation  c]  results  in  a great  simplification  feee 

reference  50  for  a more  complete  description)  of  eqns.  (15),  (15) 

(16) ,  albeit  still  not  amenable  to  analytic  solution.  In  particular, 
eqn.  (l6)  reduces  to  a Laplacian  operator  acting  on  Y with  source  terms 
on  the  right  hand  side  which  involve  the  background  ionospheric  per- 
turbed Pedersen  conductivity.  This  type  of  potential  equation  is  much 
easier  to  solve,  numerically,  than  eqn.  (l6) . Essentially,  the  barium 
cloud  motion  and  ionospheric  image  motion  evolve  on  a time  scale 

L B^/(c  E^e^) , where  L is  the  linear  dimension  of  the  cloud  perpendicular 
to  the  magnetic  field. 
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Figures  8 and  9 show  the  results  of  such  calculations.  At  t «=  0, 

b 

a barium  cloud  of  the  formZ/^  a exp  - [x^+y^J/L^  is  used  and  then  the 
small  cloud  forms  of  ecpis  (15),  (15)  (I6)  are  carried  forward  in 

time.  The  barium  cloud  drives  Images  of  order  relative  to  the  un- 

disturbed ionosphere  (see  Fig,  8).  At  long  times  the  relatively  high 
barium  cloud  ion  density  contours  (equivalent  to  Pedersen  contours 
since  v>/0  is  assumed  constant  in  each  layer),  in  the  limit  of  zero 


-1 


diffusion,  T ->  0,  form  thin  sheets  making  a definite  angle,  9 = tan 
(l^^/n^),  with  the  X cloud  velocity.  This  angular  dependence  is 
derived  from  a basic  understanding  of  the  Lagrangian  velocity  involved 
in  the  continuity  eqn  of  the  small  cloud  and  the  computer  simulations 
bear  out  this  dependence.  Figure  9 shows  the  rotational  dependence, 
in  the  plane  perpendicular  to  on  The  study  of  small  clouds 

has  given  us  great  insight  into  understanding  large  cloud  phenomenology. 
Figures  8 and  9 show  barium  cloud  morphology,  i.e,,  the  dynamical 
equilibrium  motion  (there  are  no  perturbations  put  into  the  cloud) 
which  depends  on  background  ionospheric  parameters,  in  the  small  cloud 
case,  and  image  cloud  formation.  Our  work  on  small  clouds  is  continuing 
in  an  effort  to  gain  even  more  understanding.  For  example,  we  are  now 
including  an  altitude  dependent  neutral  wind  into  the  problem  and  seeing 
the  effects  of  size  and  direction  of  such  winds  on  barium  cloud  and 
image  cloud  morphology. 

Our  large  cloud  studies  ( ^ ~ l)  use  the  full  set  of  eqns. 

P P 

(15),  (I5),  and  (16)  and  the  results  are  more  readily  interpretable 
in  light  of  our  small  cloud  studies  (see  reference  51  for  a more 
complete  description).  The  same  gaussian  form  at  t = 0 is  used  for  the 
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large  barium  cloud.  Figure  10  depicts  large  cloud  morphology  as  a 

/ b iv 

function  of  conductivity  ratio  ( S /Z'  / for  fixed  v,/0,  and  . 

'PP  bb  li 

The  case  S /7'_  =0.1  resembles  the  results  from  the  small  cloud 
^ ^ b i 

equations.  We  see  that  as  S /£  is  increased  the  cloud  rotates  more 

P P 

into  the  E^X  direction.  Indeed  as  this  ratio  gets  larger  more 
imaging  is  produced,  but  the  reaction  of  the  images  back  bn  the  motion 
of  the  barium  cloud  is  reduced.  The  results  in  the  last  row  of  Fig.  10 
are  for  an  uncoupled  cloud,  i.e,,  it  is  following  a barium  cloud  with- 
out background  ionosphere  coupling.  The  ^ = 4 case  for  the 

P P 

barivim  cloud  contours  more  resembles  the  uncoupled  cloud  case.  In  all 
cases  we  see  cloud  elongation  in  the  E X B direction  and  backside 

— O *“0 

(side  anti-parallel  to  ^ X B^) steepening.  Other  studies  show,  in 
agreement  with  small  cloud  results,  that  as  is  reduced  (to  0.1) 

the  clouds  rotate  even  more  into  the  X B^  direction  with  an  attendant 
reduction  in  strength  of  the  image  cloud. 

We  have  put  an  initial  sinusoidal  perturbation  of  about  five  wave- 
lenghts  in  the  y direction  on  top  of  the  gausslan  clouds  depicted  in 
rows  2 and  5 of  Fig.  10  to  investigate  strlation  growth  (see  reference 
31).  However,  in  terms  of  strlation  development  a one  dimensional 
cloud  is  much  more  illuminating  and  we  present  the  time  development  of 
such  a system  in  Fig.  11.  The  initial  (t  = 0)  density  Pedersen  con- 
ductivity is  given  as 


[ 


1 -3  . 10 

1 + e sin  — 

64 


p]  -0-01 


(18) 


where  iP  is  64  and  so  the  zeroth  order  cloud  is  gausslan  in  the  x 
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direction  and  infinite  in  the  y direction  (the  term  0.01  represents  a 

constant  ionospheric  background  at  the  cloud  level)  and  on  top  of  this 

is  the  initial  perturbation  e sin  (10  try/ 6k)  , i.e.,  five  wavelengths 

across  the  system.  In  the  first  row  of  Fig, 11  we  see  the  early  time 

development  of  cloud,  background  ionosphere  and  induced  potential.  As 

time  progresses  we  see  the  backside  steepening  and  backside  growth  of 

striations  and  finger-like  formation  in  the  barium  cloud  which  is  in 

2 5 

agreement  with  experimental  observations  * . What  is  of  interest  is 
the  growth  of  image  striations  in  the  background  ionosphere.  This  has 
never  been  reported  experimentally  as  observers  have  not  “looked**  for 
it.  Also  eqns  (15),  (15)  and  (l6)  in  the  ID  approximation  given  by 
(18)  are  amenable  to  a linear  stability  analysis  with  a resultant  alge- 
braic dispersion  relation,  something  that  cannot  be  performed  for  the 
2D  cloud.  Figure  12  shows  the  agreement  between  the  solution  of  such  a 
dispersion  relation  yielding  the  growth  rate  and  the  actual  numerical 
simulation  depicted  in  Fig,  11,  The  breaking  of  the  dots  (simulation) 
from  the  straight  line  (theoretical  growth)  indicates  the  onset  of  the 
nonlinear  regime. 

We  are  continuing  our  studies  on  large  clouds  for  both  ID  and  2D 
barium  clouds.  We  have  gained  a great  deal  of  insight  into  cloud  mor- 
phology and  image  formation  by  studying  its  dependence  on  ^ and 

are  currently  looking  at  the  dependence  on  We  still 

have  to  determine  onset  times  and  scale  sizes  for  striation  formation. 
Just  as  our  small  cloud  model  helped  us  to  understand  large  clouds,  we 
expect  ID  studies  to  help  us  understand  the  more  realistic  2D  clouds. 
Indeed,  the  ID  cloud  may  be  more  amenable  to  discovering  the  quasi- 
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final  saturation  or  pinching  stage.  We  also  have  to  look  at  neutral 
wind  effects  in  large  clouds  (as  we  are  doing  in  small  clouds) . We 
have  to  Include  E region  recombination  (i.e.,  chemistry)  and  its  effects 
on  image  formation. 

Most  researchers  have  focused  on  the  behavior  of  the  barium  cloud 
itself.  However,  the  major  result  of  our  two-layered  model  results  is 
the  importance  of  image  clouds,  i*e,,  coupling  of  the  cloud  to  the 
ionosphere,  under  a host  of  ionospheric  conditions.  Indeed  Images 
themselves  can  produce  anomalous  communication  and  radar  results  even 
if  our  understanding  of  the  main  cloud  is  complete.  Signals  may  have 
to  go  through  several  regions  of  turbulence  as  depicted  in  Figs.  8,  10, 
and  11  under  the  caption  cloud  and  background  ionosphere,  and  the 
regions  over  which  such  perturbations  occur  can  be  vast.  We  have  found 
that  images  always  seem  to  be  produced  (to  be  sure  their  degree  depends 
on  certain  parameters) . Although  the  basic  development  of  barium 
cloud  motion  and  striation  formation  in  these  clouds  from  our  large 
cloud  results  grossly  agrees  with  experiment  (as  does  limited  ex- 
periments on  small  cloud  agree  with  our  small  cloud  results^^),  no 
data  is  available  on  image  formation.  The  results  of  our  work  on 
coupled  plasma  cloud- ionosphere  phenomena  has  encouraged  SRI  to  look 
for  image  effects  in  a forthcoming  series  of  baritm  cloud  releases 
at  auroral  latitudes.  We  look  forward  with  great  anticipation  to 
these  experiments. 
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V,  Summary 


We  have  presented  In  Section  IV  A and  B two  classes  of  problems 

under  the  topic  of  ionospheric  irregularities.  Our  work  on  midlatitude 

nonlinear  spread  F is  proceeding  and  preliminary  results  indicate  a power 

law  for  the  power  spectrum  which  is  in  agreement  with  the  recent  in  situ 
21 

measurements  • We  wish  to  reiterate  that  the  theoretical  and  compu- 
tational programs  are  at  a stage  where  we  can  generate  detailed  infor- 
mation on  the  causes  and  structure  of  ionospheric  irregularities,  as 
manifested  in  sections  IV  A and  B,  Essentially  as  a result  of  our 
expertise  and  experience  gained  in  the  high  altitude  nuclear  problem  we 
have  been  able  to  make  significant  contributions  on  a time  scale  of 
about  one  year.  The  problems  remaining  to  be  done  can  be  accomplished 
over  a reasonable  time  period,  i.e.,  we  are  applying  existing  tech- 
nologies to  gain  a basic  understanding  of  the  causes  and  formation  of 
ionospheric  irregularities.  The  results  we  obtain  provide  a basis 
within  which  to  interpret  experimental  data, suggest  new  experiments, 
and  acts  as  a vehicle  for  suggesting  systems  design.  Understanding  the 
results  in  terms  of  the  physical  processes  involved  leads  to  a true 
predictive  capability,  i.e.,  we  can  know  what,  when  and  where  to  expect 
irregularities.  We  can  use  our  theoretical  and  computational  modeling 
in  areas  where  data  is  sparse  or  non-existent,  after  of  course  it  has 
proven  valid  in  areas  where  data  exists.  In  the  end,  this  predictive 
capability  is  very  much  needed  for  useful  systems  design. 
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The  major  emphasis  of  this  report  has  been  on  particular  types  of 

ionospheric  irregularities.  However,  in  this  concluding  section,  we 

would  like  to  enumerate  our  recent  accomplishments  in  ionospheric  model 

development  and  list  future  possibilities  (specific  future  areas  of 

study  in  ionospheric  irregularities  have  been  given  in  sections  IV 

A and  B) . The  accomplishments  are:  (1)  development  of  midlatitude 

52 

ID  ionosphere  model  ; (2)  study  of  the  role  of  resonantly  scattered 
radiation  in  maintaining  the  nighttime  E region^^;  (5)  study  of  the 
response  of  the  F region  to  EUV  flares;  (4)  proper  treatment  of  the 
role  of  oxygen  recombination  in  the  tropical  nightglow^^;  (5)  ex- 
planation of  the  non-thermal  features  of  the  auroral  plasma  due  to 

55  56 

precipitating  electrons  ; (6)  development  of  an  ionospheric  model 
for  currents,  electric  field,  and  plasma  density  in  an  auroral  arc^"^; 

(7)  development  of  the  numerical  capability  to  treat  the  coupling  be- 
tween ion  clouds  and  the  background  ionosphere,  including  imaging 
effects^^*^^;  (8)  determination  of  the  nonlinear  state  of  plasma  tur- 
bulence for  type  II  equatorial  electro jet  irregularities^^;  (9)  deter- 
mination of  the  importance  of  parallel  propagation  effects  on  the  type 

£-0 

I electrojet  instability^  ; and  (10)  development  of  a theoretical  model 
which  suggests  the  presence  of  intense  ion  heating  in  the  topside 
ionosphere  due  to  electrostatic  ion  cyclotron  turbulence^^ . 

Future  possibilities  are:  (1)  predictive  two  and  three  dimensional 

global  ionospheric  models;  (2)  predictive  two  and  three  dimensional 
global,  dynamic  models  of  the  upper  neutral  atmosphere;  (5)  three 
dimensional  models  of  the  wind  systems  in  the  polar  ionosphere;  (4) 
theoretical  and  numerical  models  of  the  ionosphere-magnetosphere 
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coupling  phenomena;  (5)  models  of  the  ionospheric  effects  of  substorms 
and  magnetic  storms;  (6)  determination,  classification  and  effects  of 
the  most  important  instabilities  which  cause  ionospheric  scintillation; 
and  (7)  generation  of  libraries  of  model  ionospheres  for  propagation 
studies  (systems  users). 

We  look  back  to  our  results  with  a feeling  of  accomplishment  and 
to  our  future  studies  with  a sense  of  excitement  and  anticipated 
accomplishment. 
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Fig.  1 — Simplified  scintillation  geometry.  A signal^ traveling  in  the  z direction  is  incident  on  a layer  of 
ionospheric  irregularities  (represented  by  ellipses).  The  scintillation  power  spectrum  in  the  observing 
plane  is  related  to  the  power  spectrum  of  the  irregularities  in  the  irregularity  layer. 
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110  km 


Fig.  2 — Model  of  daytime  equatorial  electrojet  geometry  and  observing  backscatter  radar.  The  atmo- 
spheric tidal  electric  field  E creates  a vertical  polarization  electric  field,  which  drives  the  relative 
electron-ion  drift  velocity,  V^,  by  the  Hall  effect.  This  creates  an  eastward  electrojet  current.  This 
region  of  high  conductivity  is  located  in  the  E region  ionosphere  usually  between  100  and  110  km. 

The  earth's  ambient  magnetic  field  is  taken  constant  and  in  the  x direction.  The  vertical  equilibrium 
density  gradient,  VN^,  has  a scale  length  6-10  km.  The  radar  observes  such  an  electrojet  by  pure 
backscatter  where  the  radar  wavelength  is  twice  the  wavelength  of  the  irregularity  it  views.  The  radar 
at  Jicamarca,  Peru  observes  irregularities  whose  wavelengths  are  in  the  1-10  meter  regime. 


33 


0.3 


Fig.  3 ~ Linear  stability  analysis  of  Eqs.  (2)-(6)  which  describe  the  equatorial  electrojet,  7 vs  wav^ength. 
Positive  7 indicates  growth.  For  electrojet  parameters,  Vq  ■ 100  m/sec,  ^ 4 X 10^  scc”l,  ■ 2.6  X 
10^  sec*~l  and  ® 5 X 10®  sec”!  only  waves  greater  than  28  meters  are  unstable. 
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Fig.  4 — Initial  electrojet  type  II  perturbation  at  t *.0  used  in  numerical  simulation.  These  are  the  iso- 
density contours.  The  grid  used  in  producing  these  results  was  50  X 60  mesh  poinU  denoted  by  tick 
marks  with  a 1.5  meter  mesh  spacing  in  both  directions  (y,  z),  which  models  the  electrojet  geometry  in 
Fig.  2.  Periodic  boundary  conditions  were  employed. 


34 


5 0.0 


A 


z 


m 


\ 


1 . 0 

1.0  S 0 . 0 

Y > 

Fig.  5 — Density  contours  of  type  II  electrojet  at  t 2.7  sec.  The  contours  in  Fig.  4 have  evolved  to 
this  two  dimensional  turbulent  state  via  the  electrojet  Elqs.  (2)-(6)-  The  contour  spacing  is  2.5%  of 
ambient  density  and  the  plus  and  minus  signs  represent  105  and  95%,  respectively.  This  is  the  type  of 
irregular  structure  that  propagating  signals  would  pass  through. 
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Fig.  6 — Type  II  electrojet  power  spectrum  at  t = *13  sec.  This  is  the  direction  averaged  power  spectrum 
of  the  density  fluctuations  as  depicted  in  Fig.  5.  A power  law  of  6 fairly  accurately  describes  the 
short  wavelengths.  The  fall  off  below  about  5 meters  is  due  to  grid  size  resolution  capabilities  in  the 
present  study. 
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I- ELECTRON  FLOW 

^ BETWEEN  LAYERS 


WITHIN  A LAYER 


Fig.  7 — Two  layered  model  for  plasma  cloud,  images  and  striations.  The  earth’s  ambient  magnetic  field, 
Bq*  is  ii^  the  z direction  and  taken  constant.  The  ambient  electric  field,  transverse  to  B[q.  The 

plasma  cloud  is  in  layer  1 and  image  clouds  (positive  and  negative  ion  enhancements)  in  the  background 
ionosphere  are  induced  in  layer  2.  The  model  allows  ion  flow  within  each  layer  and  electron  flow  be- 
tween each  layer  along  magnetic  field  lines.  The  model  is  two  dimensional  in  the  x-y  plane.  layer 

is  described  by  ite  own  parameters  such  as  ion-neutral  collision  frequencies,  Pedersen  conductivities,  etc., 
and  these  are  taken  to  be  constant  within  each  layer.  Typically  layer  1 is  in  the  F region  at  about  200- 
250  km  and  layer  2 is  in  the  E region  at  about  150  km. 
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Fig.  8 — Nonlinear  evolution  of  barium  cloud  Pedersen  conductivity  (Spi)  and  ionospheric  image  Pedersen 
conductivity  2*,^  for  a small  cloud: 
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2«/2:«>  = ei/2  (p,/n.)  . [1  + (^j/J2j)2]-i  (Sb/2»)  = e 2^^. 

where  2^  is  the  zeroth  order  background  ionospheric  conductivity  assumed  to  be  constant.  These  solu- 
tions are  with  no  diffusion,  T = 0,  t'j/flj  = 0.5,  the  ambient  electric  field  ^ in  the  x direction,  Bq  out  of 
the  paper  so  that  the  X velocity  frame  is  downward  (the  frame  in  which  these  results  are  depicted). 
A grid  of  128  X 128  was  employed,  with  L * 8 mesh  spaces  and  periodic  boundary  conditions.  The  inner- 
most contours  for  2^*  are  the  largest  and  for  the  images  2*  j the  innermost  contours  have  the  largest 
absolute  values  (solia  lines  depict  positive  enhancements  ana  dashes  negative  enhancements,  i.e.,  depletions, 
of  background  ionospheric  ions).  Note  that  the  scales  have  been  changes  in  the  third  row.  Dimensionless 
time  is  given  in  the  upper  right. 
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Fig.  9 — Late-time  development  of  barium  cloud  Pedersen  conductivities  for  different  values  of 
(0.0,  0.5  and  1.5).  All  other  parameters  are  the  same  as  in  Fig.  8.  As  decreases  the  cloud  lines 

up  more  with  X Bq  direction  (negative  y),  obeying  the  nonlinear  prediction  6 * tan“^  (^i/^)»  where 
6 is  the  angle  between  the  elongated  axis  of  the  cloud  and  X Bq  direction. 
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Fig.  10  — Nonlinear  dynamical  motion  of  large  2D  barium  clouds  coupled  to  the  background  ionosphere. 
This  represents  the  late  time  behavior  of  the  height-integrated  Pedersen  conductivities  for  the  barium 
cloud  and  the  induced  ionospheric  images  (left  and  center  columns,  respectively)  and  the  corresponding 
electrostatic  potential  'P  (right  column).  From  top  to  bottom  the  results  are  for  2^/Zj,  * 0.1,  1.0,  4.0 
and  oo  (uncoupled  case),  i.e.,  (X  exp  — [x2  + y2]/L2  at  t = 0 and  the  ratios  represent  the  peak 

barium  conductivity  to  unperturbed  background  ionosphere  conductivity,  at  times  t = 800  (20),  250 
(6.25),  150  (3.75)  and  350  (8.75)  sec,  respectively,  where  the  numbers  in  parentheses  are  time  in  units 
of  LB/cEq.  The  collision  frequency  parameters  are  = 0.5  and  = 0.005.  In  these  cases  ■ 

5 mV/m  and  is  in  the  negative  y direction,  Bq  » 0.5  gauss  in  the  positive  z direction  and  these  figures 
are  in  the  X B^  drift  frame,  which  is  in  the  negative  x direction.  L is  4 mesh  spaces  (4  km)  and  the 
mesh  is  32  X 32,  so  that  each  tick  mark  is  8 km.  From  the  interior  of  the  cloud  the  barium  contours 
are  e~l,  e~3,  e~^,  e"’*^  times  the  maximum  value  of  the  unperturbed  height  integrated  barium  Pedersen 
For  the  = 0.1  case  the  ionospheric  contours  are  at  1.0  ± 0.1,  ± 0.2,  ± 0.3,  ± 


conductivity, 


0.4,  ± 0.5  times  (proceeding  inward),  whereas,  for  the  Z^/Z*  * 1.0  and  4.0  cases  these  contours  are 
1.0  ± 0.1,  ± 0.3,  ± 0.5  ± 0.7,  ± 0.9  times  Z*  . The  regions  of  reWive  enhancement  and  depletion  of  the 
background  ionospheric  ions  are  marked  witn  plus  and  minus  signs,  respectively.  The  sign  of  is  indi- 
cated by  plus  and  minus  signs  in  regions  of  extrema.  These  cases  are  representative  of  F region  barium 
(plasma)  clouds  coupled  to  a background  E region  ionosphere. 
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Fig.  11  — Time  evolution  of  an  initially  ID  barium  cloud  with  an  initial  perturbation  (striations)  coupled 
to  the  background  ionosphere.  The  initial  barium  cloud  Pedersen  conductivity  is  given  by  Zp/2^  * exp 
[— x^/L^]  [l  + e“"^  sin  (10  7T  y/64)]  + 0.01  (i.e.,  2^/2^  = 1.0).  The  0.01  term  represents  a background 
F region  and  is  not  to  be  confused  with  the  background  ionospheric  E region  (where  imaging  is  occurring) 
depicted  in  the  figures.  The  geometry  and  magnitude  of  ^ and  ^ are  the  same  as  in  Fig.  10,  as  is  the 
format  for  the  figures.  In  this  case  the  collision  frequency  parameters  are  * 0.1  and  ® 0.005. 

The  times  depicted  are  t = 49.5  (0.62),  643.5  (8.04)  and  1046  (13.08)  sec,  respectively.  The  numbers  in 
parentheses  are  times  in  units  of  L Bq/cE^,.  Here  L is  8 mesh  spaces  (8  km)  and  the  mesh  is  64  X 64. 
The  barium  contours  are  now  e~i,  plus  the  background  (0.01)  times  the  maximum  value 

of  the  height  integrated  Pedersen  conductivity  of  the  cloud.  The  image  contours  are  the  same  as  the  2D 
2p/2jj  = 1.0  case.  Note  the  striation  growth  on  the  backside  (steepened  side),  i.e.,  side  away  from  direc- 
tion of  JE^  X Bo  motion.  Also  at  late  time  note  the  finger  formation  on  this  backside  (nonlinear  state). 
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Fig.  12  — Pedersen  conductivity  in  striation  region  of  cloud  versus  time  for  ID  simulation  depicted  in 
Fig.  11.  The  dots  indicate  the  values  from  the  simulation  and  the  solid  line  depicts  the  values  predicted 
by  linear  theory,  i.e.,  a linear  perturbation  analysis  of  the  full  set  of  Eqs.  (13),  (15),  and  (16).  The 
agreement  out  to  early  times  (^200  sec)  is  quite  good  and  then  the  system  takes  on  a nonlinear  character. 
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